function f = alobj(x)
    N_inequ = 4;
    global r_al pena;
    h_inequ = 0;
    g = constrains(x);
    for i = 1:N_inequ
        h_inequ = h_inequ+(0.5/pena)*(max(0,(r_al(i)+pena*g(i))).^2-r_al(i).^2);
    end
    f = obj(x)+h_inequ;
end
 %计算当前增广拉格朗日函数值
